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ABSTRACT 

* 

After  a  brief  survey  of  the  characteristics  of  a  heat  bath  and  its  role 
In  relaxation  phenomena  leading  to  the  familiar  exponential  decay.  It  is 
argued  that  the  nonexponential  form  found  commonly  in  many  condensed  matter 
systems  indicates  that  the  energy  spectrum  of  the  heat  bath  plays  a  crucial 
part  in  these  phenomena.  In  equilibrium  statistical  mechanics,  the  mean 
energy  of  a  heat  bath  determines  the  temperature  of  a  system  placed  in  contact 
with  it.  We  show  that  the  relaxation  of  a  system  placed  in  contact  with  this 
heath  bath  is  determined  by  the  distribution  of  the  energy  level  spacings  for 
level  spacings  small  as  compared  to  the  mean  spacing.  After  presenting 
arguments  in  favor  of  a  linear  behavior  of  this  distribution,  we  show,  in  a 
somewhat  heuristic  way,  that  the  resulting  relaxation  function  has  a 
nonexponential  form. 
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I.  INTRODUCTION 

Experimental  data  on  relaxation  phenomena  In  diverse  areas  of  condensed 
matter  physics  are  quite  generally  found  to  exhibit  slower  than  exponential 
decqy  for  long  times  in  the  form  {see  ref.  1.  and  references  therein) 

exp  C-a(t/Ts  )bL  a>o,  o<b<l  (1) 

where  tsis  a  characteristic  time  in  the  system.  Traditionally2*3  the  residual 
part  of  relaxing  quantities  are  discussed  either  in  terms  of  a  pure 
exponential  decay 

exp[-t/t  ]  (2) 

where  T  is  the  "relaxation  time", or  in  terms  of  a  superposition  of  such  terms 
with  a  distribution  of  T.  ts*  t  in  general  but  there  is  a  relation  between 
them  (see  Ref.  1.).  (References  4,5,6  are  a  representative  set  selected  at 
random  out  of  a  large  number  of  papers  on  this  subject.  See  Ref.  1.  for  more 
citations).  In  this  framework,  t  is  obtained  from  a  calculation  of  the  time- 
independent  transition  rate  given  by  the  Golden  Rule.  The  physical  picture  of 
relaxation  here  is  that  the  system  which  is  relaxing  and  which  is  described  by 
a  Hamiltonian,  H$,  is  in  contact  with  a  heat  bath  which  is  a  much  larger 
system  described  by  a  Hamiltonian,  Hg,  and  which  is  not  affected  by  the 
interaction  (Hg$)  with  the  relaxing  system.  This  interaction  is  supposed  to 
be  "weak"  and  leads  to  relaxation.  The  equilibrium  properties  of  the  given 
system  in  contact  with  the  heat  bath  are  determined  by  the  temperature  of  the 
bath,  which  in  turn  is  just  the  mean  energy  associated  with  the  bath.7*8 

An  alternative  scheme  to  derive  (2)  is  the  master  equation  approach  (see 
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for  example  Ref.  8).  Here  again  the  time-independent  transition  rate  is 
employed  in  setting  up  the  master  eqatlon.  The  physical  model  for  relaxation 
Is  however  still  the  same  -  system,  heat  bath,  and  their  mutual  (weak) 
Interaction.  An  elementary  discussion  of  the  exponential  decay  in  stochastic 
processes  and  in  quantum  mechanics  may  be  found  in  Merzbacher's  book,9  and  an 
extensive  discussion  of  the  role  of  time  independent  transition  rates  Is  found 
in  Tolman7. 

Now,  since  exponential  decay  Is  not  observed,  and  the  decay  law  given  by 
(1)  Is  more  a  rule  than  an  exception,  it  is  natural  to  seek  an  explanation  for 
this  behavior  by  examining  In  more  detail  the  origin  of  the  time-independent 
transition  rate.  This  shift  of  emphasis  from  time  Independent  transition  rate 
(T1TR )  to  time  dependent  transition  rate  (TDTR)  In  order  to  arrive  at  (1)  has 
been  emphasized  recently  by  Teltler  et.al.10  from  phenomenological 
considerations  of  rate  equations,  and  from  general  considerations  based  on  the 
Paley-Wlener  theorem  by  Ngal  et.al.11  It  was  also  recognized  early  in  the 
development  of  the  TDTR  that  using  the  Golden  Rule1  required  a  linear 
dependence  of  the  level  spacing  density  for  low  spaclngs  in  order  to  obtain 
the  long  time  behavior  of  the  form  given  by  Equ.  (1).  The  purpose  of  the 
present  paper  is  to  suggest  that  the  heat  bath  be  described  by  any  chaotic 
quantum  system.  The  reasons  for  this  suggestion  are  given  In  Sec.  II.  In 
Sec.  Ill,  we  set  up  the  calculatlonal  scheme  for  computing  TDTR  and  obtain  the 
required  time  dependence  In  terms  of  the  slope  and  cut  off  of  the  linear  law 
for  the  level  spacing  distribution  at  small  spaclngs  of  the  heat  bath.  In 
Sec.  IV  a  brief  summary  of  the  results  obtained  Is  given. 
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II.  SPECTRAL  CHARACTERISTICS  OF  THE  HEAT  BATH 

We  are  primarily  Interested  In  the  long  time  relaxation  of  any  physical 
property  of  a  system,  for  example  a  dielectric  or  a  mechanical  property.  What 
is  Involved  In  a  relaxation  process  Is  a  readjustment  of  the  system  in  such  a 
way  that  there  Is  no  transport  phenomenon  accompanying  It.  In  general 
terms,  the  system  that  Is  undergoing  relaxation  Is  described  by  a  Hamiltonian 
Hj,  and  it  Is  supposed  to  be  in  contact  with  a  much  larger  system  described  by 
a  bath  (or  reservoir)  Hamiltonian  Hg.  The  interaction  between  the  two,  whose 
Hamiltonian  is  Hgg,  Is  assumed  to  be  weak  and  Is  supposed  to  Induce  the 
process  of  relaxation  In  the  physical  quantity  of  Interest  In  the  system.  It 
Is  Important  to  realize  that  the  bath  system  Is  large  compared  to  the  system 
that  is  under  investigation  so  that  while  the  bath  is  not  affected  by  the 
Interaction  with  the  system.  Its  effect  on  the  system  is  paramount. 

The  precise  nature  of  the  heat  bath  is  left  unspecified  except  for 
stating  that  the  system  In  contact  with  It  acquires  its  temperature  in  the 
equilibrium  situation.  In  the  conventional  approach,  the  details  of  the  bath 
hamlltonlan  are  not  Important  even  though  It  Is  recognized  to  have  an  almost 
continuous  energy  spectrum,  by  virtue  of  its  enormous  size.^  Also,  in  the 
final  analysis,  the  bath  variables  do  not  occur  In  the  description  of  the 
system  so  that  one  averages  over  these  variables  In  computing  properties  of 
the  system.  Since  only  the  temperature  of  the  bath  enters  the  picture  the 
only  relevant  entity  appearing  In  this  picture  is  the  mean  energy  of  the  bath 
system,  which  Is  kept  fixed,  thus  determining  the  temperature.  It  is  clear 
from  this  description  of  the  heat  bath  that  a  detailed  knowledge  of  Hg  Is  not 
required,  except  for  Its  temperature  and  the  obvious  observation  that  Hg  has 
:;st  continuous  spectrum  bounded  from  below  with  a  finite  mean.  One  of 


the  approaches  to  the  theory  of  relaxation  Is  to  construct  the  density  matrix 
associated  with  the  system  plus  bath  In  the  presence  of  Hg$  and  Integrate  out 
the  bath  variables  by  tracing  over  all  the  states  of  the  bath,  leaving  behind 
a  residual  density  matrix  for  the  system,  which  can  be  used  to  compute  any 
property,  P(t),  of  the  system,  by  calculating  the  appropriate  average  of  the 
operator  representing  P (t )  over  the  residual  density  matrix.  Feynman  and 
Vernon13  have  given  a  formal  path  Integral  representation  for  this  density 
matrix  and  Fano1^»15  has  given  an  expression  in  the  Interaction 
representation. 

It  may  not  be  out  of  place  here  to  point  out  that  there  are 
circumstances  when  the  nature  of  the  heat  bath  Is  known  purely  from  the 
physical  consideration  of  the  energy  or  time  domain  one  is  examining.  For 
certain  electronic  properties,  phonons  (the  motion  of  the  crystal  lattice)  are 
the  relevant  heat  bath  system  and  so  for  excitations  Involving  frequencies  of 
the  order  of  1013  Hz  and  above,  the  phonon  excitations  determine  the 
relaxation  properties  upto  times  of  order  10"13  sec  and  here  the  relaxation 
rate  Is  essentially  exponential.  When  one  waits  for  longer  times,  say  10-10 
sec,  the  phonon  excitations  are  not  relevant  any  longer  and  the  relevant  bath 
system  must  be  something  else.  It  must  be  pointed  out  however,  that  the 
phonon  system  is  Itself  imbedded  In  the  new  bath  system  which  we  are  proposing 
so  that  there  Is  a  common  temperature  for  all  the  entities  making  up  the 
system.  Thus  we  may  picture  a  hierarchy  of  heat  baths,  each  Imbedded  in  the 
other,  so  that  they  all  have  a  common  temperature  and  each  Is  a  relevant  heat 
bath  In  the  appropriate  time  regime.  This  nesting  of  heat  baths  leads  to 
different  time  dependences  of  the  relaxing  entity  in  different  time  domains. 
What  we  are  Interested  In  for  the  present  work  is  the  relatively  long  time 
domain  such  as  10"1®  sec  and  lower,  where  the  usual  known  excitations  become 


irrelevant  and  a  new  mechanism  is  called  for.  In  the  present  paper  we  are 
concerned  only  with  this  regime.  It  must  be  stressed  that  the  theoretical 
formalism  given  by  Fano1**16  is  applicable  quite  generally  to  all  these 
situations. 

In  the  traditional  description^  no  mention  of  the  nature  of  the  heat  bath 
energy  spectrum  is  made  except  for  its  being  continuous.  Thus  we  describe  the 
system  undergoing  relaxation  by  means  of  the  Hamiltonian 

H  =  H$  +  Hjj  +  Hb$  (3) 

where  HB,  HB,  and  HBS  are  the  hamiltonian  for  the  system,  bath,  and  their 
mutual  interaction.  One  may  then  compute  the  density  matrix  of  the  entire 
system,  given  that  at  time  t*0  the  system  and  the  bath  are  not  interacting, 
and  have  been  prepared  such  that  the  system  Is  in  some  preassigned  state  and 
the  bath  is  in  thermal  equilibrium. 

As  a  model  for  the  heath  bath  we  shall  adopt  any  large  quantum  system 
whose  classical  motion  Is  Irregular.  For  such  systems  Berry  has  shown  in  a 
series  of  elegant  papers16"20  that  the  quantum  levels  are  fairly  regularly 
distributed  and  that  the  probability  density  P(S)  for  the  spacings  between 
neighboring  levels  has  the  asymptotic  form 

P(S)«  aS  (0<S«5)  (4) 

provided  S  Is  small  as  compared  to  the  average  spacing  5  .  For  larger  values 
of  S  the  spacings  distribution  P(S)  goes  through  a  maximum  to  decay  to  zero  at 
large  values  of  S.  These  details  depend  on  the  precise  choice  of  the  quantum 
system,  but  the  linear  behavior  of  the  spacings  distribution  appears  to  be 
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universal,  i.e.  It  holds  for  a  “generic"  chaotic  quantum  system. 

It  is  remarkable  that  the  same  linear  behavior  (4)  of  the  spaclngs 
distribution  Is  found  if  the  heat  bath  Is  described  by  means  of  a  random 
matrix  hamlltonlan.  The  hamiltonlan  for  the  heat  bath  is  very  complex  so  that 
we  may  replace  It  by  a  statistical  description.  For  the  determination  of 
equilibrium  properties  of  a  system  in  contact  with  the  heat  bath,  only  the 
average  of  the  bath  hamiltonlan  is  needed.  Hence  one  can  try  to  use  a 
"Gaussian  Orthogonal  Ensemble"  (GOE)  for  the  random  matrix  for  describing  the 
heat  bath  because  we  take  the  bath  system  as  being  time  reversal  invariant. 
Since  only  the  mean  value  of  the  bath  Hamiltonian  is  required  for  equilibrium 
properties,  we  use  a  "canonical  ensemble"  in  setting  up  its  density  matrix. 

We  now  observe  that  the  GOE  has  known  average  spectral  properties,  which  we 
employ  in  our  analysis  of  the  Golden  Rule  in  determining  the  TDTR  induced  by 
the  bath  in  the  system.  For  a  description  of  the  philosophy  and  the  theory  of 
random  matrix  hamiltonlans  one  may  refer  to  Porter's  collection  of  papers  and 
his  clear  introductory  summary21,  Mehta's  book22,  and  a  more  recent  review  by 
Brody  et.  al.23. 

It  is  remarkable  that  both  models  for  the  heat  bath  which  we  have 
considered  in  this  section,  (a)  an  irregular  quantum  system  and  (b)  a  random 
matrix  hamiltonlan,  lead  to  the  same  linear  behavior  of  the  spacings 
distribution.  We  feel  that  (a)  is  the  physically  correct  model  for  a  heat 
bath  and  that  the  random  matrix  is  just  a  convenient  way  to  simulate  an 
irregular  quantum  system.  It  may  not  be  out  of  place  here  to  conjecture  that 
the  random  matrix  hamiltonlan  may  Indeed  be  a  very  accurate  model  for  a 
generic  Irregular  quantum  system  If  one  coarse  grains  the  energy  spectrum.  A 
hint  of  this  equivalence  may  be  found  when  one  compares  the  average  density  of 
states  for  the  quantum  version  of  Sinai's  billiard  19  with  the  middle 
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part  of  the  semi-circular  law  appropriate  for  the  density  of  states  of  a 
rarw'om  matrix. 21-23 

III.  CALCULATION  OF  THE  TIME  DEPENDENT  TRANSITION  RATE 


The  calculation  proceeds  in  three  steps:  (i)  set  up  a  rate  equation  for 
the  physical  quantity  that  is  undergoing  relaxation  in  terms  of  a  TDTR;  (ii) 
compute  the  TDTR;  and  (iii)  solve  the  rate  equation  once  the  TDTR  is 
determined.  Our  goal  is  to  examine  the  long  time  limit  of  the  time  dependence 
and  so  there  is  much  simplification  that  can  be  made  right  from  the  start.  A 
formal  justification  of  step  (i)  is  given  elsewhere24.  We  are  concerned  here 
mainly  with  step  (ii)  and  calculate  the  transition  rate  using  the  Golden  Rule 
with  proper  attention  paid  to  the  types  of  interactions  that  could  be 
suggested  for  HBS.  Having  done  this,  we  then  invoice  the  cumulant  expansion 
technique14  or  equivalently  the  linked  cluster  scheme  to  calculate  the  TDTR 
essentially  to  all  orders  In  Hgj.  This  is  in  the  same  spirit  as  in  the  binary 
correlation  approximation. 

We  are  thus  led  to  consider  the  survival  probability  of  a  state  of  the 
system  when  it  Interacts  with  the  bath.  Let  the  heat  bath  have  states  |b>  and 
the  relaxing  system  two  representative  states  |1>  and  |2>.  At  t*o,  the 
combined  system  (S  +  B)  is  In  a  state  |s,b>  where  |b>  is  some  state  with  its 
energy  In  the  range  E  -AcE^xE  +  A  .  For  a  “good"  heat  bath,  in  general,  one 
does  not  know  the  state  |b>  apart  from  the  fact  that  it  has  an  energy  Eb  in 
some  energy  region  (E-  A*  E+  A).  The  standard  Golden  Rule  result**  for  the 
probability  that  the  total  system  Is  in  the  state  ( s * b ' >  at  time  t  given  the 
initial  state  is  |sb>  is 


10 


|Csv(t)|2  *  2  |VS|jtS'5'|^  (1-cos  “sb.s'b't)  w  ^sb.s'b* 

r 

where  |i(  wsb>s«bi  =  Es*+  Eb<  -  Es  -  Eb,  with  Es,  the  system  energy.  The  matrix 
element  of  the  system-bath  coupling  Hamiltonian  Hgj  between  the  states  |sb> 
and  | s  * b ‘ >  is  denoted  by  Vsb# s  * b *  - 

In  order  to  obtain  the  total  transition  probability  QSjS*(t)  for  the 
system  to  go  from  state  |s>  tz  state  | s * >  irrespective  of  the  state  of  the 
heat  bath  we  sum  over  |b‘>  and  average  over  |b>.  This  gives  for  the 
probability  to  find  the  system  in  state  J s ' >  at  time  t 


s,s 


.u)  -  V 

111 


?,b'lVsb,  s'b'l2  (1"cos  “sb.s'b'1^  “sb.s'b”  (5) 


where  B  denotes  the  number  of  heat  bath  states  in  (E-a,E+a).  In  order  to 
express  the  relaxation  function  Qs>s*(t)  in  terms  of  the  spectral 
characteristics  of  the  heat  bath  one  now  proceeds,  in  a  somewhat  heuristic 
fashion,  as  follows. 

Firstly,  as  the  system-bath  coupling  will  depend  only  on  the  global 
properties  of  the  heat  bath  the  matrix  element  Vsb>sibi,  will  be  practically 
independent  of  the  choice  of  |b>  and  | b * >  in  the  energy  range  (E-&,  E+A). 
This  enables  us  to  put 


Vsb,s'b' 


(6) 


and  to  bring  the  constant  VSts'  outside  the  summation  sign  in  (5). 

Secondly,  we  note  that  the  resulting  sum  over  b'  will  be  independent  of 
the  choice  of  b.  Hence  the  average  over  b  is  trivial  and  (5)  can  be  written 
as  the  sum 
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Qs,s,(t)=  'jj?  ^s.s1'2?  (1"cos“sb,s,b,t)  “sb.sV 


(7) 


In  a  relaxation  process,  the  transition  occurs  from  a  state  of  the  system 
to  another  state  of  the  system  which  is  essentially  degenerate  with  it  and  so 
N“sb,sV  a  Eb'  “  Eb*  The  sum  (7)  can  now  be  written  as  an  integral 


Qs,s>(t>  ■  4 


I*/ 


(1-  COS  et/K)  e’2  P(e)  de. 


(8) 


where  p(c)  de  denotes  the  average  number  of  heat-bath  quantum  states  with 
energies  in  (E^+  e,  E^+  e+de)given  a  heat-bath  level  at  E^.  For  elarge  as 
compared  to  the  average  level  spacing  6one  has  P(e)~  6  On  the  other 
hand,  for  e  <  4  we  can  use  the  result  (4)  for  the  distribution  of  level 
spacings  of  an  irregular  quantum  system  and  put  P(e)~  «  e. 

Thirdly,  we  use  the  following  qualitative  form  for  p  suggested  by  the 
preceeding  remarks 

p(e)  *  a  e  0<e<EQ,  (9a) 

p  (e)  =  6”1  cq  <  e  <  •,  (9b) 


e  « 
o 


(9c) 


The  calculation  of  the  Integral  (8)  is  now  straightforward,  and  leads  to  the 


asymptotic  behavior 

Q  (t)“  4  <*|  V  .|2  *n4)  + .  (t»««W),  (10a) 

t  *  a  4  M  exp(-l-T),  (10b) 
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The  dots  In  (10a)  denote  terms  that  vanish  for  t  and  the  constant  y=0.5772 
is  Euler's  constant.  Of  course*  for  small  values  of  t  the  quantity  Qs>$i(t) 

O 

will  be  proportional  to  tS 

The  calculation  outlined  above  amounts  to  lowest  order  perturbation 
theory.  It  was  shown  by  Fano1*  that  when  one  proceeds  in  a  rigorous  fashion 
one  obtains  a  cumulant  expansion  for  the  transition  probability,  which 
essentially  leads  to  an  exponentiation  of  the  Golden  Rule  result  (10) 

QSfS.(t)  ~  (t/T)b.  (t»t),  (11a) 

b  *  4  ot|V$  s.|2.  (lib) 

With  a  time  dependent  transition  rate  of  this  form  the  solution  of  a  rate 
equation  will  lead  to  terms  of  the  form  (1),  with  0<b<l,  for  otherwise,  with 
b>l,  one  has  faster  than  exponential  decay. 

IV.  CONCLUDING  REMARKS 


We  have  shown,  in  a  somewhat  heuristic  way,  how  the  fine-grained  spectral 
characteristics  of  the  heath  bath  determine  the  form  of  the  relaxation 
function.  It  Is  remarkable  that  the  only  quantities  which  enter  are  the 
average  spacing  6,  the  average  system-bath  matrix  element  VS(S«,  and  the 
slope  oof  the  spacing  distribution  at  small  spaclngs.  The  linear  behavior  of 
this  distribution  at  small  spaclngs  Is  a  generic  feature  of  Irregular  quantum 
systems16"20,  hence  an  Irregular  (chaotic)  quantum  system  Is  a  universal  model 
-  heath  bath. 
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Our  considerations  also  show  that  one  can  use  the  Gaussian  orthogonal 
ensemble  of  random  matrices  to  simulate  a  chaotic  quantum  system,  and  hence  to 
model  a  heath  bath.  This  might  also  explain  the  success  of  random  matrix 
theory  in  nuclear  physics  and  other  branches  of  physics. 

The  main  motivation  of  the  work  reported  In  this  paper  was  to  make  the 
concept  of  the  heat  bath,  which  up  till  now  Inhabited  the  literature  on 
statistical  physics  as  an  almost  featureless  entity,  more  specific  and  to 
determine  which  of  its  properties  enter  into  the  physics  of  relaxation. 
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